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Abstract  -  Radar  offers  unique  advantages  over  other  sensors 
in  the  human  detection  problem,  such  as  remote  operation  during 
virtually  all  weather  and  lighting  conditions.  Many  radar-based 
human  detection  systems  today  employ  Fourier  analysis,  such  as 
spectrograms.  However,  spectrograms  perform  poorly  in  high 
clutter  environments.  Also,  an  inherent  SNR  loss  is  caused  by  the 
implicit  assumption  of  linear  target  phase.  In  this  paper,  human 
modeling  is  used  to  derive  a  more  accurate  non-linear 
approximation  to  the  true  non-linear  target  phase  and  the 
liklihood  ratio  is  optimized  over  unknown  parameters  to  enhance 
detection  performance.  Performance  is  compared  both 
analytically  and  through  MATLAB  simulations. 


I.  Introduction 

Much  of  the  current  research  in  human  detection  is 
focused  on  the  application  of  Fourier  transform  based 
techniques,  such  as  spectrogram  analysis.  In  2002, 
Geisheimer  [1]  experimentally  showed  that  a  measured 
human  spectrogram  could  be  obtained  by  summing  the 
spectrograms  obtain  from  the  reflections  off  each 
individual  body  part.  This  result  was  theoretically 
confirmed  by  Van  Dorp  [2],  who  based  his  calculations 
on  a  kinematic  walking  model  developed  by  Thalmann 
[3].  Comparisons  between  the  human  spectrogram  and 
the  spectrograms  of  other  animals,  such  as  dogs,  showed 
that  certain  characteristics  were  unique  to  just  humans 
and  that  these  characteristics  could  be  used  in  human 
detection  and  identification.  For  example,  Otero  [4] 
proposed  techniques  based  on  features  extracted  from 
the  spectrogram  to  differentiate  humans  from  dogs,  or 
even  women  from  men.  Greneker  [5]  designed  and 
tested  a  suicide  bomber  detection  system  based  on 
variations  in  the  spectrogram  caused  by  the  presence  of 
a  bomb. 

However,  all  of  these  systems  were  used  in  close 
proximity  or  in  situations  where  the  clutter  was  minimal. 
In  practice,  spectrograms  are  almost  completely 
obscured  in  high  clutter  environments,  rendering 
impossible  the  accurate  extraction  of  any  information 
[6].  Furthermore,  application  of  the  FFT  implicitly 
assumes  that  the  target  phase  history  is  linear,  an 
assumption  not  generally  true  for  non-linear  phase 
human  targets.  The  resulting  phase  mismatch  results  in 
an  inherent  SNR  loss  due  to  the  detector,  decreasing 


even  further  the  chances  of  detecting  slow-moving 
human  targets. 

In  this  paper,  Thalmann' s  kinematic  model  is  used  as 
a  basis  for  deriving  a  more  accurate  non-linear 
approximation  to  the  true  target  phase.  Maximum 
liklihood  estimates  (MLE)  of  unknown  geometry  and 
target  parameters  are  obtained  to  maximize  the 
likelihood  ratio  and  resulting  matched  filter  output.  The 
performance  of  the  FFT -based  matched  filter  is 
compared  to  that  of  the  proposed  optimized  linear  and 
nonlinear  phase  detectors,  as  well  as  to  the  ideal 
"clairvoyant"  detector,  representing  the  best 
performance  attainable  with  complete  knowledge  of  the 
target. 


1 1 .  Signal  modelling 

In  general  the  received  radar  signal  is  comprised  of 
the  sum  of  noise  plus  a  time-delayed  and  phase-shifted 
version  of  the  transmitted  chirp  signal.  For  the  purposes 
of  this  paper,  the  noise  is  modeled  as  being  complex 
Gaussian.  Then  the  target  return  may  be  expressed  as 

sr(n,t)=atrect\^-  0) 


where  at  is  the  amplitude  as  given  by  the  radar  range 
equation;  *  =  2 R/c ;  R  is  the  range  from  the  antenna  to 

target;  c  is  the  speed  of  light;  x  is  the  pulse  width;  y  is 
the  chirp  slope;  fc  is  the  transmitted  center  frequency; 
and  n  is  the  pulse  number. 

If  the  noise  is  small  compared  to  the  signal  energy, 
the  peak  of  the  pulse  compression  output  occurs  at  the 
range  bin  in  which  the  target  is  present: 


,M= 


-- 


-j^Rln] 


R2[n] 


(2) 


where  A  is  an  unknown  amplitude  parameter  and  the 
slow-time  dependence  on  pulse  number,  n,  is 
emphasized  with  brackets.  Typically  this  peak  is 
indiscernable  in  noise,  so  that  the  detection  test  must  be 
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applied  separately  over  each  range  bin.  However,  in  the 
case  of  human  targets,  the  center  of  the  range  bin,  rb,  is 
not  a  sufficiently  accurate  approximation  to  the  true 
target  location.  Most  human  target  motion  will 
generally  remain  within  one  range  bin,  excepting  the 
case  of  very  long  dwell  times. 

A  better  approximation  of  target  range  may  be 
obtained  by  assuming  that  the  motion  is  linear  along  a 
constant  angle,  0,  relative  to  the  initial  antenna-target 
vector  (see  Fig.  1).  First  define  it  as  the  vector  between 
the  aperture  and  initial  target  location;  r„  as  the  vector 
between  the  aperture  and  target  location  at  the  nth  pulse; 
and  h  as  the  vector  describing  human  motion.  Then, 


In  this  work,  the  human  body  has  been  divided  into 
twelve  parts:  head,  torso,  two  upper  arms,  two  lower 
arms,  two  thighs,  two  lower  legs  and  two  feet.  Each 
body  part  is  modeled  as  a  point  target  with  a  cylindrical 
radar  cross  section  (RCS),  excepting  the  head  which  is 
modeled  with  a  spherical  RCS.  The  total  radar  return  is 
given  by  the  sum  of  the  returns  from  each  point  target. 

Each  kinematic  equation  depends  on  two  variables: 
1)  RV,  the  ratio  of  velocity  (v)  to  thigh  height  (HT), 
and  2)  t%,  a  time  index  taken  relative  to  the  beginning 
of  a  step.  Mathematically, 

RV  =  —  (6) 

HT 


rn\Z  =  Irl | 2  +  H 2  -  2|rl ||h| cos  d 


(3) 


Since  |h| « |rx|  and  ^  \  +  */  for  small  x, 


lrnHrlHhlcos0 


(4) 


Final  Target 
Location 


Figure  1 .  Target  Geometry 


A.  Linear  Phase  Approximation 

If  we  model  Ihl  as  simply  a  constant  velocity,  v,  times 
the  pulse  repetition  interval,  T,  times  the  pulse  number, 
n,  then  subsituting  (4)  into  (2)  results  in  a  linear  phase 
model 

,  (5) 

rb 

where  the  range  term  of  the  amplitude  has  been  modeled 
as  rb  since  the  effect  of  perturbations  in  the  amplitude  is 
negligible;  and  liqhr  an  unknown  parameter. 

B.  Non-Linear  Phase  Approximation 

The  Thalmann  kinematic  model  for  human  walking 
consists  of  a  set  of  equations  that  describe  the  time 
varying  position  and  change  in  angle  of  key  joints  and 
limbs.  Combined,  they  can  be  used  to  calculate  the 
time-varying  position  of  any  point  in  the  human  body. 


and 


t%=J?L^+t  (7) 

where  t0  is  a  constant  indicating  the  point  within  the 
stepping  cycle  that  the  first  transmitted  pulse  reflects 
from  the  target.  The  velocity  profile  of  each  body  part 
can  be  clearly  seen  in  the  resulting  human  spectrogram. 
The  strongest  return  results  from  the  low-amplitude 
sinusoidal  motion  of  the  torso,  while  the  largest 
Doppler  amplitude  is  given  by  the  trajectory  of  the  feet. 
A  typical  measured  spectrogram  can  be  visually 
matched  to  a  sample  simulated  spectrogram,  as  shown 
in  Fig.  2. 


(a)  (b) 

Figure  2.  Human  spectrograms:  (a)  Measured  (from  [4]) 

(b)  Simulated. 

Because  of  the  strong  signal  strength  of  the  return 
from  the  torso,  only  the  torso  equations  of  motion  will 
be  used  in  calculating  Ihl: 

|h|2  =  OSy  +  osl  +  osjg  (8) 


where 


OSv  =  0m5RV[-\  +  sm{\m%-0.1Jt)} 
OSL  =  Al  sin(2®%  -  0.2#) 


J  -0.032 
0. 1287?  V  2  +0.128 


for  0.5  <  if  V  <  2.3 
for  RV  <  0.5 


and 
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OSFB  =  Aa  sm(4m%  +  2 </)a ) 

f  -0.021  for0.5<RF<2.3 

a  _{-0.084W2_0.084W  for  RV  <0.5 
(f>a=  -0.127  +  0.731>/W 

The  torso  model  can  be  further  simplified  by 
neglecting  the  case  when  RV  <  0.5,  as  most  human 
targets  will  not  fall  in  that  range  and  by  then  averaging 
the  trajectories  for  0.6  <  RV  <  2.3,  the  maximum  RV 
for  which  Thalmann’s  model  is  defined.  The  individual 
torso  trajectories  and  resulting  average  trajectory  is 
shown  in  Fig.  3. 


Individual  and  Average  Torso  Trajectories 


Figure  3.  Individual  and  Averaged  Torso  Trajectories. 


A  sinusoidal  curve  fit  to  the  average  torso  trajectory 
may  be  then  found  to  be 

|h|  =  h[n] «  0.0108cosf^^V^V«  +  ^j  + 0.037  ^ 


estimate,  ,  which  will  change  depending  on  whether 

we  apply  an  FFT,  or  a  linear  or  nonlinear  approximation 
to  Ihl. 

Assuming  that  we’ve  obtained  a  suitable  model  for 
the  expected  target  return,  a  likelihood  ratio  test  may  be 
formed  to  determine  the  detector  decision  rule. 


zN  R]  II 


p{x;Hi)  = 


_  1  „-(x-y . )”  Rj‘  (x-y . ) 


r lR, 


Decide  H.if  p[x:  Re{v»  R  / 

0) 


(11) 


where  /  =  ^Re)yj)lR|1yp|  Q~l(PFA)  and  pfa  is  the 
desired  probability  of  false  alarm. 

A  FFT -Based  Detector 

Taking  the  FFT  across-slow  time  for  each  range  bin 
yields  the  following  estimate  to  the  true  target  signal: 


Kif, 


ypJM' 


p  V  peak  > 

Nil 


-j2jtfPe*kTn 


(12) 


B.  Optimized  Linear  Phase  ( OLP)  Detector 

Inspection  of  (5)  reveals  three  unknown  parameters 
over  which  the  likelihood  ratio  must  be  maximized:  A, 
r,  and  Vr  =  vcos0-  If  we  express  yp=As,  the  an  MLE 
estimate  of  A  may  be  found  as 

\  _  Re{s"x}  (13) 


where  RV  and  §  are  unknown  parameters.  This 
expression  may  be  substituted  into  (2)  and  (4)  to  obtain 

r  A  ~j~^L{r~^n\  COS  0) 

ypL»J*-2«  c 

(10) 


III .  Detector  design 


The  unknown  parameters  in  the  phase  may  be 
estimated  from  the  tangent  of  the  phase,  found  by  taking 
the  ratio  of  the  imaginary  and  real  parts  of  the  signal. 
This  operation  will  convert  the  complex  Gaussian  noise 
w  into  Cauchy  noise  with  distribution 


p(x-,n,r)= 


(14) 


The  detector  must  make  a  decision  between  two 
hypotheses: 

H0 :  x  =  w 
Hi  :  x  =  y  p  +  w 

where  w  is  Gaussian  noise  with  covariance  matrix  Ri 
and  yp=yp[n]  is  the  true  target  signal.  As  it  is  not 
possible  to  know  apriori  the  interference  covariance 
matrix,  instead  an  estimate  is  used,  Rr  =  £[xH0x||0]  • 
Similarly,  the  true  target  signal  is  approximated  with  an 


Here,  p  is  a  location  parameter  for  the  Cauchy 
distribution  that  is  yp  under  Hi  and  zero  under  H0.  The 
parameter  y  is  an  unknown  constant  that  drops  out 
during  maximization. 

The  MLE  estimates  of  r  and  vr  can  be  shown  to  be 
found  by  solving  the  following  equation,  where  the 
derivative  is  taken  with  respect  to  r  or  vr,  depending  on 
which  estimate  is  being  computed: 


I 

n= 1 


x[n]  - 


Im{y„W}Ylm{y„[n]} 


Re{y„[n]}l  Re{y„[n]} 


(15) 
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The  resulting  set  of  three  equations  with  three 
unknowns  is  highly  nonlinear,  making  impossible  a 
closed  form  expression  for  the  MLE  estimates. 
However,  these  equations  may  be  numerically  solved  in 
MATLAB  to  derive  a  numerical  solution  for  the 
estimates.  When  the  MLE  estimates  A,?,andvr  are 
substituted  into  (5)  and  used  in  (1 1),  the  OPT  detector  is 
obtained. 

C.  Optimized  Non-Linear  Phase  ( ONLP)  Detector 

The  same  procedure  as  outlined  in  the  previous  section  is 
employed  to  derive  the  MLE  estimates  for  the  optimized  non¬ 
linear  phase  detector  based  on  the  nonlinear  motion 
approximation  expressed  in  (9).  With  the  Thalmann-based 
expression  for  h,  (15)  is  numerically  solved  in  MATLAB  to 
obtain  MLE  estimates  for  the  unknown  parameters  A,  r,  VR, 
<|>,  and  0.  Substitution  of  these  estimates  into  (10)  and  (11) 
yield  the  ONLP  detector. 

D.  Clairvoyant  Detector 

Naturally  the  best  possible  performance  is  obtained  if  we 
have  sufficient  information  to  identically  know  the  true  signal 
embedded  in  noise,  yp.  Thus,  the  optimal  clairvoyant  detector 
is  obtained  by  using  yp  -  yp  in  (11). 


IV.  Performance 

Detector  performance  is  compared  both  analytically  and 
via  simulations  in  MATLAB.  The  expected  receiver 
operating  characteristic  (ROC)  curves  may  be  mathematically 
written  as 

PD  =  Q^-1(PFA)-^L{ffR^\}\  (16) 

Simulations  were  conducted  for  a  human  target  walking  at 
2  m/s  approximately  7  km  away  from  a  radar  antenna  with  a 
transmitter  power  of  1.8  kW,  center  frequency  of  1  GHz,  and 
a  dwell  time  of  100  ms. 

For  this  example,  the  estimates  of  the  phase  of  yp  are 
compared  in  Fig.  4.  Notice  that  the  nonlinear  phase 
approximation  closely  follows  the  true  target  phase. 

Theoretical  and  simulated  PD  versus  PFA  and  PD  versus 
SNR  curves  are  shown  in  Figs.  5,  6,  7  and  8.  As  expected, 
the  nonlinear  phase  detector  outperforms  all  other  detectors, 
performing  just  under  the  ideal  clairvoyant  detector.  At  an 
SNR  of  -20  dB,  the  ONLP  detector  has  a  PD  that  is  roughly 
five  times  higher  than  that  of  the  FFT -based  detector.  The 
SNR  loss  for  the  ONLP  detector  is  4  dB  less  than  that  of  the 
FFT-based  detector  for  most  PDs.  The  OLP  detector  also 
performs  much  better  than  the  FFT-based  detector,  although 
still  falling  short  of  the  ONLP  detector. 


Figure  4.  Approximated  and  True  Signal  Phase  Histories 


SNR  =  -20  dB 


Figure  5.  Theoretical  ROC  Curves  for  All  Detectors 
at  a  SNR  of  -20  dB 


SNR  =  -20  dB 


Figure  6.  Simulated  ROC  Curves  for  All  Detectors 
at  an  SNR  of  -20  dB 
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V.  Conclusions 


PFA  =  0.01 


Figure  7 .  Theoretical  Probability  of  Detection  versus  SNR 
at  a  1%  False  Alarm  Rate  for  All  Detectors 


SNR  (dB) 

Figure  8.  Simulated  Probability  of  Detection  versus  SNR 
at  a  1%  False  Alarm  Rate  for  All  Detectors 


A  new  human  target  detector  design  is  presented 
based  on  more  accurate  modeling  of  the  true  target 
phase  history  and  optimization  of  the  likelihood  ratio 
test  is  presented.  Theoretical  analysis  and  MATLAB 
simulations  of  performance  show  tliat  the  new  detector 
achieves  substantially  higher  detection  probabilities 
over  a  range  of  SNRs  relative  to  the  currently  obiqutous 
FFT -based  detector. 
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